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Abstract 

A new line list including positions and absolute intensities (in the form of Einstein A values and oscillator strengths) has been 
produced for the OH ground X 2 n state rovibrational (Meinel system) and pure rotational transitions. All possible transitions are 
included with v' and v" up to 13, and .1 up to between 9.5 and 59.5, depending on the band. An updated fit to determine molecular 
constants has been performed, which includes some new rotational data and a simultaneous fitting of all molecular constants. The 
absolute line intensities are based on a new dipole moment function, which is a combination of two high level ab initio calculations. 
The calculations show good agreement with an experimental v=l lifetime, experimental /a values, and Av=2 line intensity ratios 
from an observed spectrum. To achieve this good agreement, an alteration in the method of converting matrix elements from 
Hund’s case (b) to (a) was made. Partitions sums have been calculated using the new energy levels, for the temperature range 
5-6000 K, which extends the previously available (in HITRAN) 70-3000 K range. The resulting absolute intensities have been 
used to calculate O abundances in the Sun, Arcturus, and two red giants in the Galactic open and globular clusters M67 and M71. 
Literature data based mainly on [O I] lines are available for the Sun and Arcturus, and excellent agreement is found. 
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1. Introduction 

The OH radical is very important in atmospheric chem¬ 
istry due to its high reactivity with volatile organic com¬ 
pounds mm, and it is the major oxidizing species in the lower 
atmosphere 0. OH is also produced in the upper atmosphere in 
excited vibrational levels, and near infrared emission to lower 
levels is the main cause of the nighttime airglow of the up- 
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per atmosphere @0. This airglow is often an unwanted fea¬ 
ture in astronomical observations 0, but has sometimes been 
used for wavelength calibration 0 and atmospheric tempera¬ 
ture retrievals 0. OH is also present in stars 0, interstellar 
clouds 0, extra-terrestrial atmospheres MM, and often in 
relatively large concentrations m in flames 0! El US- 

The transitions of interest in this work are in the Meinel sys¬ 
tem, which are the rovibrational transitions within the X 2 n 
ground state. These have been used to determine the oxy¬ 
gen abundance in the Sun lfl6l and other stars, for example by 
Melendez and Barbuy flTl and Smith et al. Ifl8l . 

There have been several ab initio dipole moment functions 
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Figure 1: Calculated and experimental dipole moment functions for the OH 
X 2 n state. 


(DMFs) calculated for the OH ground state (T9J [20] j21] [22 :, 
1221241. Accurate experimental dipole moments were obtained 
for the v=0, 1, and 2 levels by Peterson et al. l25l . and these 
were used by Nelson et al. 127] l26ll in combination with their 
own experimental relative line intensities to calculate a DMF 
at internuclear distances between 0.70 and 1.76 A. Another 
experimental DMF was obtained by Turnbull and Lowe 


also using the Peterson et al. values and their own experimen¬ 
tal intensities. The various DMFs have all disagreed to some 
extent l26l (also see Figure [I]), including around the equilib¬ 
rium internuclear distance ( r e ), where the intensities are most 
sensitive to the DMF. 


A set of transition probabilities was calculated by Mies |f29l , 
and more recently Goldman et al. l30l produced a set of inten¬ 
sities, which is now the most widely used for abundance calcu¬ 
lations, for example by Melendez and Barbuy ITT) . This list by 
|Goldman et al.| is currently in the HITRAN Oil database, and 
more extensively in the HITEMP | [32l database . The intensities 
are based on the DMF of Nelson et al. l26l between 0.70 and 
1.76 A, and that of Langhoff et al. If22l at other internuclear dis¬ 
tances. Another line list including ab initio Einstein ,4 /- / values 
was calculated recently 02 l24l . however the accuracy can be 
improved, and it did not cover all of the available vibrational 
levels (up to v=9 as opposed to v=13). This DMF is included 
in the comparisons in Figure [T] 

Bernath and Colin 132 performed a fit to obtain molecular 
constants up to v=13 in 2009, using the best data available at the 
time. Lines from the extensive study of Melen et al. 1341 were 
included, along with the rovibrational measurements of Abrams 
et al. (151 (22 bands, Av=l-3, up to v=10), Nizkorodov et al. 
1351 (11-9 band), and Sappey and Copeland 1361 (12-8 band), 
and B 2 S + -X 2 n system lines from Copeland et al. 1371 (v=7- 
13 in the X 2 n state and v=0-1 in the B 2 X + state). They also 


assigned new rotational lines in v=4 in an ATMOS (1994) 1381 
solar spectrum. 

The fit of Melen et al. 1341 included rotational lines in v=0- 
3, from both their assignments of an ATMOS (1985) 1391 solar 
spectrum and laboratory measurements, and rovibrational labo¬ 
ratory data from Maillard et al. m. Hyperfine rotational data 
in v=0-3 from several other studies l40ll44ll42l 142l44l 1431 l46l 
|42|42|42|32|3I]|52l|32|5ll52|36l|37l|38l were also included 

in the fit by Melen et al. l34l . Despite hyperfine structure not 
being included in the results of Bernath and Colin 132 . they 
took some of this information into account by adding pseudo¬ 
transitions to the fit, between e and / parity levels and various J 
levels (“A-doubling data”), based on the term values calculated 
by Melen et al. l34l . For a further description of the sources 
of the included transitions, the reader is referred to Table 1 of 
Bernath and Colin 132 and Figure 4 of Melen et al. 1341 . 

Due to limitations in the fitting program used by Bernath and 
Colin 1331 . four separate fits were required to calculate all of the 
molecular constants. In addition, more rotational data in v=0- 
2 has recently been reported l59l . which has a higher reported 
accuracy than the previously available lines. A fit in which this 
new data was included and all parameters were floated simulta¬ 
neously was performed in this work, and is described in Section 

m 

In this paper, the currently available HITRAN intensities are 
compared with experimental data (a v= 1 lifetime l60l . Av=2 H- 
W ratios (this work), and /z v values {25]). Based on these com¬ 
parisons, it was concluded that the calculation of a new set of 
line intensities would be beneficial, and Sections[2to|2describe 
their calculation and validation. Partition function calculations 
are described in Section [7] and Section [8] uses the newly cal¬ 
culated line intensities to calculate the O abundance in the Sun 
and three other stars, to demonstrate their use in astronomical 
calculations. 

2. Molecular Constant Fit 

A new fit of molecular constants was performed, which was 
based mainly on that of Bernath and Colin 8331 . Line weight¬ 
ings from the previous fit were retained, all parameters were 
floated simultaneously, and lines were replaced by new pure 
rotational data from Martin-Drumel et al. lf59l where possible, 
which was the only new data included. The A-doubling data in¬ 
cluded by Bernath and Colin 1331 was also included here, and as 
with their fit, no hyperfine structure constants were calculated. 

The fit of Bernath and Colin 13311 was very extensive in terms 
of the number of molecular constants obtained. A number of 
higher order constants that had unreasonable fitting uncertain¬ 
ties were fixed at values based on extrapolation and careful in¬ 
spection of the energy level patterns. 

The same molecular constants as used by Bernath and Colin 
132 were used here, and those that were fixed previously were 
also fixed. Floating these constants was attempted, but their 
resulting uncertainties were too large. As all parameters were 
now being floated, the uncertainties in general increased, and 
some of the higher order constants that were successfully de¬ 
termined in the fit of Bernath and Colin 132 were fixed at their 
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previous values. Compared to the work of Bernath and Colin 
El. the small amount of new data and simultaneous determi¬ 
nation of all constants resulted in only a small improvement in 
the fit. Specifically, the total unweighted average line position 
error changed from 0.02649 to 0.02637 cm -1 , and the weighted 
average error improved by a factor of 1.32. 

See Table[I]for the final molecular constants. 
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Table 1: Molecular constants 3 for the OH X 2 n state (in cm *). 


V 

T 


A 

B 

D x 10 3 H x 10 7 

Lx 10“ 

0 

0 


-139.050877(41) 

18.53487308(79) 1.9089352(64) 1.42658(14) 

-1.4707(11) 

l 

3570.35244(72) 


-139.32031(62) 

17.82391990(96) 1.870447(10) 1.37737(33) 

-1.5700(35) 

2 

6975.09355(69) 


-139.58940(54) 

17.1224805(18) 1.835828(12) 1.31875(33) 

-1.6991(41) 

3 

10216.1478(12) 


-139.84424(86) 

16.427914(26) 1.80599(16) 1.2447(37) 

-1.824(32) 

4 

13294.6144(14) 


-140.0832(17) 

15.737018(34) 1.78267(21) 1.1630(37) 

-2.124(19) 

5 

16210.5997(18) 


-140.2893(19) 

15.045630(63) 1.76640(58) 0.9956(143) 

-1.416] 

6 

18962.9921(19) 


-140.4396(23) 

14.348665(68) 1.76159(67) 0.8425(176) 

-1.75] 

7 

21549.1777(19) 


-140.5041(19) 

13.639345(70) 1.77073(73) 0.6013(204) 

-1.74] 

8 

23964.6390(21) 


-140.4296(16) 

12.90888(12) 1.8002(18) 0.3214(773) 

-3.4] 

9 

26202.4452(19) 


-140.1409(16) 

12.145349(64) 1.85782(47) [-0.111] 

-5.8] 

10 

28252.5370(31) 


-139.5064(17) 

11.33258(17) 1.9571(20) [-0.717] [- 

-12.5] 

11 

30100.8094(29) 


-138.3238(17) 

10.45025(33) 2.1383(85) [-1.5] 


12 

31728.2544(35) 


-136.3625(19) 

9.45859(52) 2.349(13) [-2.5] 


13 

33109.6605(30) 


-133.0318(23) 

8.32780(25) 2.7906(46) [-4.0] 


V 

Mx 10 15 

N x 10 20 OxlO 23 

y X 10 

7D x 10 5 

y« x 10 9 y L x 10 13 

qx 10 2 

0 

0.8842(25) 

[6.2] 

[-3.96] 

1.19251(14) 

2.357(13) 

-3.021(218) 3.47(103) 

-3.86908(17) 

l 

1.258(11) 

[-8.80] 

M.O] 

1.13749(54) 

2.315(18) 

-2.467(282) [3.5] 

-3.69357(18) 

2 

1.721(14) 

[-32.6] 

M-0] 

1.07690(48) 

2.192(15) 

-1.523(228) [3.5] 

-3.51941(29) 

3 

2.007(92) 

[-66.0] 

M-0] 

1.0250(32) 

2.353(159) 

-1.96(111) [3.5] 

-3.3393(21) 

4 

[2.5] 

[-100] 

H.O] 

0.9664(42) 

2.473(274) 

-3.14(252) [3.5] 

-3.1590(33) 

5 

[2.8] 

[-100] 


0.9035(46) 

2.591(266) 

[-6.4] 

-2.9719(40) 

6 

[3.1] 

[-100] 


0.8390(56) 

3.191(361) 

[-17.1] 

-2.7800(45) 

7 

[3.4] 

[-100] 


0.7593(51) 

3.240(362) 

[-17.7] 

-2.5781(49) 

8 

[3.7] 

[-100] 


0.6744(63) 

4.353(576) 

[-30.0] 

-2.3593(75) 

9 

[4.0] 

[-100] 


0.5558(64) 

5.313(648) 

[-60] 

-2.1184(82) 

10 

[4.3] 

[-100] 


0.4010(85) 

[6.0] 

[-100] 

-1.861(29) 

11 




0.2981(190) 

[7.0] 


-1.564(68) 

12 




0.4404(207) 

[8.5] 


-1.284(31) 

13 




0.8677(115) 

[10.0] 


-0.7239(207) 

V 

q D X 10 5 

q H X 10 9 q L X 10 13 

q M X 10 17 

P 

Pd x 10 5 ph x 10 9 

PL x 10 12 

0 

1.4739(14) 

-2.713(29) 4.358(225) 

-3.42(26) 

0.235355(20) 

-5.216(21) 5.635(380) 

-1.050(194) 

l 

1.4451(14) 

-2.591(32) 3.961(133) 

[—2.08] b 

0.224684(30) 

-5.235(27) 7.775(529) 

[—2.88] b 

2 

1.4342(15) 

-2.786(29) 5.180(125) 

[-2.0] 

0.213190(81) 

-5.003(30) 4.386(460) 

[—3.20] b 

3 

1.403(13) 

-2.505(201) 4.748(886) 

[-2.0] 

0.20276(35) 

-5.310(233) 6.54(191) 

[-3.0] 

4 

1.400(19) 

-2.484(137) [4.79] b 

[-2.0] 

0.19120(43) 

-5.296(260) [3.80] b 

[-3.0] 

5 

1.388(21) 

-2.5] 



0.17937(59) 

-5.701(438) [3.5] 


6 

1.413(26) 

-2.75] b 



0.16642(59) 

-5.949(485) [2.8] 


7 

1.427(30) 

-2.69] b 



0.15223(63) 

-6.479(566) [1.8] 


8 

1.464(66) 

-3.0] 



0.13578(86) 

-7.111(957) [1.0] 


9 

1.475(79) 

-3.3] 



0.11694(88) 

-8.29(116) 


10 

1.612(394) 

-3.6] 



0.09316(231) 

-8.96(459) 


11 

[l-7] b 




0.04981(193) 

[-12.000] 


12 

[1.8] b 




0.03033(168) 

[-18.000] 


13 

[1.9] b 




-0.05720(168) 

[-27.000] 


a 

Constants in square 

Drackets 

were fixed; one standard deviation in 

the last digits is 

given in parentheses; the constants with their full 


precision as used by PGOPHER are available in the supplementary material. 



b 

Fixed here at the value from Bernath and Colin 1331, but floated in their fit. 
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3. Kitt Peak Experimental Spectrum and Herman-Wallis 

Ratios 

Relative intensities were extracted from an experimental FTS 
spectrum from Abrams et al. US. for later comparison to 
our calculated values. The spectrum was recorded at the Kitt 
Peak National Observatory, Arizona, by observing OH emis¬ 
sion from the H + O3 reaction, and it has an excellent signal- 
to-noise ratio. The intensity axis of the observed spectrum 
was calibrated using a NIST traceable tungsten ribbon lament 
lamp (Optronics Model 550C), a spectrum of which had been 
recorded using the same equipment. The Av=2 region was cho¬ 
sen for analysis as the OH lines are strongest here, it is the 
cleanest region (mostly due to the presence of water lines in 
the Av= 1 region), and the calibration curve is much flatter in 
the Av=2 region. Lines were fitted using wspectra ED, the 
areas were taken as the intensities, and Herman-Wallis (H-W) 
ratios (described below) were obtained for bands up to 9-7. 

OH contains a light H atom, and so exhibits a strong H-W ef¬ 
fect li62l , meaning that the transition intensity heavily depends 
on 7. The obvious effect on the spectrum is that the P branches 
are more intense than the R branches. When comparing cal¬ 
culated intensities to an observed spectrum, this effect can be 
exploited to validate the line parameters. The intensity ratio 
of an R branch line to a P branch line with the same upper 7 
level will vary with 7, and as the same upper (for an emission 
spectrum) 7 level is used, the effects on intensity due to energy 
level population are canceled. This intensity ratio is hereafter 
referred to as the H-W ratio. 


4. Line Intensity Calculation Method 

4.1. Summary 

For a full description of the general line intensity calculation 
procedure, see our other recent work E3ll64ll . as these calcula¬ 
tions were carried out in the same manner. 

Briefly, equilibrium constants (Table [2} were calculated by 
least-squares fitting to the usual power series expansions in v 
+ 1/2 (equations shown in the supplementary material), using 
the constants 7 V and B w from Table [l] These were employed 
in the program RKrIES to generate a potential energy curve, 
which was then entered into LEVEL li66l along with the DMF. 
LEVEL generates vibrational wavefunctions and transition ma¬ 
trix elements (MEs), but does not include electron spin. These 
MEs are therefore in terms of N and not J, and we refer to 
them as Hund’s case (b) MEs. They were converted to case 
(a) MEs which include electron spin using a transformation 
equation, which was derived in our recent CN work |63l , and 
adjusted here in Section |4.2| PGOPHER li67l sets up standard 
N 1 Hamiltonians using the molecular constants from Table [I] 
and adjusts the case (a) transition MEs to include the required 
averaging over the angles between space and body fixed axis 
systems. The Hamiltonians are then diagonalized, and the re¬ 
sulting eigenvectors are combined with the transition matrices 
to give ’’transformed” transition matrices in terms of the real 
states. The square of these MEs are the line strengths, S j'j, 


Table 2: Equilibrium molecular constants (in cm -1 ) for the OH X ’ [ I state. 
These constants with the full precision used in the calculations are available in 
the online supplementary material. 


Constant 3 

Value 

Constant 

Value 

oj e 

3738.465(19) 

B e 

18.894867(49) 

LJ e X e 

84.875(18) 

Vei 

0.72343(15) 

u e y e 

0.5409(77) 

(Xe 2 

0.007212(138) 

£O e Ze 

-0.02252(170) 

a e3 

-0.0006656(469) 

UeTJ ei 

-0.0009854(1986) 

a e4 

0.00005108(598) 

Uerie 2 

0.00003087(1155) 

a e 5 

-0.000004828(251) 

u e ile 3 

-0.000004539(264) 




a Numbers in parentheses indicate one standard deviation to the 
last significant digits of the constants. 


which can be used in the following equation to calculate Ein¬ 
stein A values ||68l : 


16n 3 v 3 S j’j 
3e 0 hc 3 (2J' + 1) 

r>3c 

= 3.136 188 94 x 1CT 7 —— J -^~. 

(27' + 1) 


( 1 ) 

( 2 ) 


We report our line intensities in the form of Einstein A values, 
and also oscillator strengths (/-values), which are calculated as 
follows: 


m e eoC 3 (27' + 1) 

(3) 

27re 2 v 2 (27+1) Aj ' j 

1 (27' + 1) 


1.499 1938— -A ri . 

(4) 

v 2 (27 + 1) 


where S j'j is in debye, and Aj-j is in seconds. 


4.2. Transformation from Hund’s case (b) to (a) Matrix Ele¬ 
ments 

For transition strengths, we require the MEs of the space 
fixed electric dipole operator, T p (ji) in spherical tensor notation. 
This is conventionally expanded lf69l in terms of a molecule 
fixed dipole moment operator, T k (/j), and a Wigner D matrix, 
D k pq {u> ), that gives the transformation between the space and 
molecule fixed axis systems: 

7 >) = X D iM Tk M (5) 

In a Hund’s case (a) basis the wavefunctions can be expressed 
as a product of a rotational part, \JMEl), and a vibronic part, 
\rjAv;SI.), where rj represents the remaining electronic and vi¬ 
brational quantum numbers. Taking matrix elements of the 
dipole operator, the standard approach yields a product of 
terms: 

^(/M'O'ID^MI/MOX^A'tA'S'ir^llTyA^S), (6) 

9 

with the first term containing integrals over the rotational wave- 
functions leading to Honl-London factors and the second term. 
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the vibronic transition moment being a band strength indepen¬ 
dent of rotation. To allow for case (a)/(b) mixing the standard 
methods will give rotational wavefunctions that are linear com¬ 
binations of the basis functions, giving slightly more compli¬ 
cated expressions. 

Given a potential energy curve, V(r), derived as described 
above, the vibrational wavefunctions, v F(r), can be calculated 
by solving the one dimensional Schrodinger equation: 

-ti 2 d 2y VJj) 

+ KM'FvW = E'i’yir), (7) 

The vibronic transition moment is then evaluated by integrating 
the DMF, /j(r) over the vibrational wavefunctions involved: 

J f*00 

y v'(r)/j(r)'¥ v (r)dr ( 8 ) 

o 

The H-W effect arises because the vibrational wavefunctions 
change with rotation which can be modelled by adding a cen¬ 
trifugal term: 


and ignoring the electron spin in finding the vibrational wave- 
functions. The centrifugal term becomes B(r)(^N(N + 1) - A 2 ) 
where N is the angular momentum excluding electron spin 
(J - N±\ here). This will give different results; the most obvi¬ 
ous change is that the effective case (a)/(b) mixing term involves 
an average value of B, rather than an operator that depends on 
r. We expect the differences to be small, given that this is any¬ 
way the approach used in calculating energy levels, and there 
are less stringent requirements on the accuracy of intensity cal¬ 
culations. Integrating the dipole moment over wavefunctions 
derived in this manner gives essentially vibronic transition mo¬ 
ments in a case (b) basis, which must be transformed to a case 
(a) basis for use by PGOPHER, as it always works in a case (a) 
basis. The required transformation has been derived and used 
in our previous papers, but it became clear as part of this work 
that a small but important correction is required, so we now 
summarise the revised derivation of the transformation. 

The transformation between case (a) and functions, defined 
in terms of Q = A + X, and case (b) functions, defined in terms 
of N, is fTP : 


2nr 


(/(/ + 1) - Q 2 ) = B{r)[j(J + 1) - O 2 ) 


(9) 


to the potential energy in the Schrodinger equation above. A 
separate vibrational wavefunction, V F V j{r) then arises for each 
value of /, and the formally vibronic transition moment ac¬ 
quires a variation with J. To reflect this we add additional labels 
to the transition moment: 

J f*00 

^jdrUr^J^dr (10) 

o 

This means that the Schrodinger equation must be solved sepa¬ 
rately for each value of J. 

There are two approaches to taking the spin into account. In 
the method described by Chackerian et al. EH, that has been 
used by Goldman et al. l30l to prepare the line list used in HI- 
TRAN 2012, a pair of coupled differential equations are set up, 
one for Q = 1/2 and one for Q = 3/2. This essentially corre¬ 
sponds to a Hund’s case (a) basis. The diagonal terms are the 
one-dimensional Schrodinger equation as above, with the addi¬ 
tion of a constant spin-orbit coupling term, ,4 AX = ±A/2. Note 
that the spin-orbit coupling constant. A, was assumed to be in¬ 


dependent of r by Chackerian et al. This is an approximation. 


though table [I] indicates a relatively small variation of A with 
v and thus r. The off-diagonal term is -B(r) \JJ(J + 1) - GO', 
analogous to the - BJ+S _ term normally used to model case 
(a)/(b) mixing (“spin uncoupling”). Solving the pair of equa¬ 
tions gives a vibrational wavefunction for each v, J, and Q, and 
the vibronic part of the transition moment is then evaluated by 
integrating the DMF over each pair of wavefunctions on the 
assumption that it is independent of fi. While the dipole mo¬ 
ment itself should be diagonal in Q (with case (a) basis func¬ 
tions), transitions off-diagonal in O will arise by virtue of the 
case (a)/(b) mixing. 

We use a simplified approach here, allowing us to use 
LeRoy’s LEVEL program to solve a single differential equation. 


I 77 A; NAS JM) = 

^(- 1 y v - s+A+£ V2A + J ( A + 2 |t/A; SX; JMD). 

( 11 ) 

This allows any given case (a) matrix element to be expressed 
in terms of case (b) MEs: 

J> 1 )V( 2 N' + 1 )( 2 A + 1 ) (£ 

q ' ' 

X (o -X -^a) W^'-’N'A'S'JM\H\tiA;NASJM ). 

( 12 ) 

The matrix element of the transition dipole moment given in 
equation (2) above in a Hund’s case (a) basis is: 

(77'A; SX'; J'M'£l'\Tp(jS)\i)A-,S?,', JMQ) = ^(-l) M '~ n ' 


x V(2/' + l)(2/+l )|_‘ T ' a 

X (rj'A'\Tq(J'£l'JQ.)\i)A). 
and in a Hund’s case (b) basis: 


/ J' 

k 

A 

( r 

k 

j \ 


q 

nj 

\-M' 

q 

m) 


(rj'A; N'S J’ M'\T k p (qi)\qA\ NS JM) = 

X (-M' \ m) (~\) N +SJrJ+k y/(2 Jr + l)(2J + 1) 

X {T n l } Z ( - 1 )W, " A ' V^' + IX^+D 

X (-A' q 


(13) 


(14) 
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Table 3: Symmetrised e transition matrix for the (2,0), R(1.5) example tran¬ 
sition level. Upper - with original transformation equation. Lower - including 
AS =£ 0 MEs. All values are in debye. 



l 2 n,i.5)0//)) 

l 2 n (0 .5)O//)) 

l 2 n, i. 5) o//)> 

-0.01157 

0 

l 2 n (0 .5)O//)> 

0 

-0.01418 


l 2 n ( i. 5) 0 //)) 

l 2 n (0 .5)O//)> 

l 2 n (1 . 5 ) 0 //» 

-0.01157 

-0.0002774 

l 2 n ( o.5)0//)> 

-0.0001114 

-0.01418 


In this expression, (tj'A'\T^(N'N)\t]A) is the vibronic tran¬ 
sition dipole moment which only depends on N by virtue 
of the H-W effect, analogous to the J and Q dependence in 
(rfA'\Tq(J'Q.'JQ.)\r]A). (q' A'\T^(N' N)\qA) is the quantity cal¬ 
culated by LeRoy’s LEVEL program. Combining the last three 
equations allows the Hund’s case (a) and Hund’s case (b) vi¬ 
bronic matrix elements to be related: 


WA'\T k n ,_ n (J'Q.'m)\r 1 A) 


\j'~ fi' 


J' k J 
-Q! q Q 


(-D 

x (IN' + V)(2N + 1) 


z<-» r 

N,N' 

J' S N' 
O' -I 


~N+n'-n+s+j+A'+k 


J 


\N' 

J' 

S \ 

IN' 

k 

N\ 

\J 

N 

k ) 

l-A' 

q 

a) 


A 'Mo 

k 


N 

-A 


(15) 


A more complete derivation is given in the supplementary ma¬ 
terial. The key difference from the previous version l64l is the 
removal of the condition that 0 = 0', and in the OH case this 
results in matrix elements with AO = ±1 in addition to the 
AO = 0 matrix elements expected in a simple Hund’s case (a) 
basis. That this is a direct consequence of the N (or J) depen¬ 
dence of the vibronic transition moment follows from the result 
(shown in the supplementary information) that if the vibronic 
transition moment is independent of N the expected result fol¬ 
lows: 


(q'A'\T k A ,_ A (N'N)\qA) = < 7 'A'|r*,_ n (y'fl'70)|i 7 A>, (16) 

and the vibronic transition moment is independent of J and has 
the selection rule AO = A A (or A2 = 0). 

The H-W effect therefore has the effect of introducing small 
matrix elements off-diagonal in O in the case (a) basis. For 
example, for the (2,0), R(1.5) transition, the transition matrix 
changes as shown in Table [3] 


4.2.1. Comparison to HITRAN Calculations 

To show that the results obtained using this adjusted trans¬ 
formation equation are reasonably equivalent to the potentially 
more accurate HITRAN style calculation involving separate vi¬ 
brational wavefunctions for each Q component, calculations 
were performed using a DMF equivalent to that used for HI¬ 
TRAN, and also the HITRAN molecular constants. The DMF 


was constructed using the Nelson et al. ||26|| DMF at short range 
(the same range as used for HITRAN), and our ACPF DMF 
(Section [5j at long range (as the long range part used for HI¬ 
TRAN does not appear to have been published). This means 
that compared to the actual DMF used for HITRAN, it is iden¬ 
tical at short range, but then diverges slightly. Hereafter, it is 
referred to as the ’’HITRAN DMF”. 

Figure [2] plots the ratio of the transition intensities obtained 
using both the revised and original transformation methods, 
and the HITRAN DMF and molecular constants. Much bet¬ 
ter agreement is seen with the use of the revised transformation 
method. As the potential used in this comparison is not exactly 
that used for HITRAN, some of the discrepancies are due to to 
this small difference. Transitions are only shown up to v' = 2 
as the potential and DMF will begin to diverge from those of 
HITRAN at higher v levels. 

The improvement in the transformation method is further 
demonstrated in Figure [3] which compares the H-W ratios ob¬ 
tained with the previous transformation method and HITRAN 
DMF, with those observed as described in Section[3] The main 
feature in this graph is the splitting of the solid (An transitions) 
and dotted ( F 22 transitions) lines of the same color. This split¬ 
ting is caused by the distribution of intensity between the F\ \ 
and F 22 transitions. The splitting pattern of the red lines clearly 
does not match that of the observed lines, but the splitting of the 
HITRAN lines matches well, although the HITRAN lines are 
displaced upwards on the graph from the observed lines (which 
is caused by inaccuracies in the DMF). With the use of the re¬ 
vised transformation equation, the red lines would match the 
blue HITRAN lines almost exactly. 
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10' 3 10' 2 10' 1 10 ° 10 1 10 2 10 3 
Einstein A (this work; s' 1 ) 


Figure 2: Ratio of Einstein A values between HITRAN and obtained from us¬ 
ing the methods described in this work but with the HITRAN 2012 molecular 
constants and DMF. a: Using the previous version of the transformation equa¬ 
tion. b: using the final version of the transformation equation. Transitions 
are shown for all bands (including pure rotational) with V up to 2, and where 
F' = F". “Observed” means that the transitions have been identified in experi¬ 
mental spectra, but does not refer to observed intensities. 


—■— HITRAN 
♦ Our calculations 
—•— Observed 


Solid lines - FI-FI 
Dotted lines - F2-F2 



Figure 3: H-W ratios in the OH X 2 n (2,0) band. The values plotted are equal 
to the R branch intensity divided by the P branch intensity for transitions that 
share an upper J and F level. The red lines are calculated using our methods 
(but the previous version of the transformation equation) and potential, and the 
HITRAN DMF and molecular constants. 


5. Dipole Moment Functions 


5.1. Final DMF and Ab Initio Calculations 

A new DMF was used in this work, which is a combination 
of two ab initio DMFs. At large internuclear distances, an aver¬ 
aged coupled-pair functional (ACPF) 1721173117411751 DMF was 
used, and at small distances, a unrestricted open shell coupled- 
cluster method with single, double, and perturbative triple exci¬ 
tations (UCCSD(T)) l76l DMF was used. 

The ACPF DMF was calculated using MOLPRO 2012.1 1771 . 
using an aug-cc-pV6Z basis set GO |79]|8()]]. This followed a 
complete active space self-consistent field (CASSCF) ffSTl [82l 
calculation using an active space of four <x, four n. and one 
6 orbital. The UCCSD(T) DMF was calculated with MOLPRO 
2010.1 ll76l . and with the aug-cc-pV6Z basis set IfTHi [79 801. 
For both DMFs, all core correlation was included, and the 
dipole moments were calculated by the finite field method. 

A comparison of these DMFs and others is shown in Figure 
[7] A full description of the mixing of the two DMFs to produce 
the final version is provided in Section 5.3 


5.2. Results From Using Only Ab Initio DMFs 

The full calculations were performed using both ab initio 
DMFs individually, and the results were analysed by compari¬ 
son to the experimental dipole moments of Peterson et al. H25I . 
an experimental lifetime obtained recently by van de Meerakker 
et al. l60l . and H-W ratios. 

Better agreement with the Peterson et al. l25l dipole moments 
was shown with the ACPF DMF than any of the other DMFs 
shown, except for that of Nelson et al. which is expected as 
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Table 4: Calculated and experimental lifetimes of the OH X 2 n, v=l, 
7=1.5, F i level, using different DMFs. The first column is the literature 
source of the DMF. The second is the lifetime reported in that paper, and 
the third is the lifetime calculated using the specified DMF, but with our 
methods, potential, and line positions. 


DMF 

Lifetime (ms) 

Reported by Calculated in 

authors this work 

van de Meerakker 

et al. (2007) (60l (expt) 

59 ±2 


van der Loo and 

Groenenboom (2008) 12311241 

56.84 

56.97 

Langhoff et al. (1989) 1221 

73.3 

66.7 

HITRAN 2012 QT) 

56.6 ± 10-20% 

56.6 

This work 
(ACPF DMF) 


65.36 

This work 
(UCCSD(T) DMF) 


59.14 

This work 

(Final DMF; UCCSD(T)+ACPF) 


59.13 


Figure 4: Calculated and experimental // v values of OH (X 2 n). The error bars 
of the experimental values for v=0 and v= 1 are not shown as they are slightly 
smaller than the size of the symbols. 


the Nelson et al. DMF is partly based on these experimental 
measurements. This is shown in Figure [4] where it can also be 
seen that the UCCSD(T) DMF does not agree as well. 

The lifetime of the v= 1 level was calculated as the weighted 
(for e and / parity levels) sum of the Einstein A values for all 
possible transitions from a single upper J level (7=1.5, Ff), 
and taking the reciprocal. Unfortunately, the result of 65.36 ms 
does not fall within the error bounds of the recent experimen¬ 
tal measurement. The UCCSD(T) DMF produces a lifetime in 
excellent agreement of 59.14 ms. Lifetimes were also calcu¬ 
lated using our methods, potential, and line positions, but using 
the other literature DMFs. All of these are shown in Table [4] 
which shows that the UCCSD(T) lifetime is the best match to 
the experimental value. 

H-W ratios calculated from the Kitt Peak spectrum as de¬ 
scribed in Section [3] were obtained for Av=2 bands up to (9,7), 
and compared to our calculated values and those of HITRAN. 
Results from the two DMFs are quite similar for the low vi¬ 
brational bands, and are closer to the observed ratios than HI¬ 
TRAN. This is shown for the 2-0 band as an example in Figure 
[6] For the higher vibrational bands, the calculated ratios show 
worse agreement with the observed ratios for the UCCSD(T) 
DMF than the ACPF DMF. 


UCCSD(T) method does not perform well at longer range (in 
this case above about 2.1 A). Therefore, it was decided to use 
the UCCSD(T) at short range, the ACPF at long range, and to 
mix them in the intermediate region. 

The best intermediate region that gave a resulting DMF with 
smooth first and second derivatives was found to be 1.5 to 2.1 
A. Also to achieve this, a value was added to the ACPF DMF. 
To ensure that the ACPF still decreased smoothly to zero, this 
value was smoothly decreased to zero with increasing range. 

Specifically, the ACPF DMF was split into three regions, 
one where the full amount was added (r < 1.5 A), one where 
the amount added was smoothly decreased to zero (1.5 A < 
r 3.0 A), and one where nothing was added (r > 3.0 A). The 
amount added in the smoothing region was equal to 

cos j —(r - 1.5) 

where a is the maximum amount added. 

In the mixing region, the final DMF is a linear combination 
of the two adjusted ab initio DMFs: 

DMF = c(ACPF) + (1 - c)(UCCSD(T)), (18) 
where the coefficient c is equal to 

^(c°s — l'5)j + lj- ( 19 ) 




5.3. Mixing of Ab Initio DMFs to Form Final DMF 

The v=l lifetime is determined mainly by the gradient of the 
DMF around R e , and as can be seen in Figure [T] this varies 
between the various DMFs that have been calculated. The 
shape of the DMF also determines the magnitude of the H- 
W effect, and as the UCCSD(T) DMF reproduces these well, 
its shape is likely to be accurate. The match to the experi¬ 
mental yU v values can be improved by subtracting a constant 
from the DMF, whilst retaining its shape. The single-reference 


A fit of the DMF was performed, in which the final intensi¬ 
ties were fitted to the Av=2 H-W ratios (from 100 upper J 
levels) and the Peterson et al. l25l /a values, and the param¬ 
eters floated were the maximum value added to the ACPF DMF 
(0.04988 D), and the value subtracted from the UCCSD(T) 
DMF (0.001534 D). The /r v values were weighted such that 
their weights were proportional to reciprocal of their reported 
uncertainties, and the sum of their weights was equal to the sum 
of the weights of all of the H-W ratio data. The H-W ratios were 
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Figure 5: Mixing of the two ab initio DMFs to form a final DMF. 


weighted using the reciprocal of their uncertainties (calculated 
as described in Section[3]i. The v=l lifetime was not used in the 
fit, as all of the floated parameters have a negligible effect on it, 
and it is already reproduced well by the DMF. The mixing of 
the two DMFs using the final parameters is shown in Figure [5] 
and the DMFs are all available in the supplementary material. 

5.4. Comparison to Nelson et al. 4261/ Fit 

Nelson et al. l26l calculated their DMF by fitting an expan¬ 
sion about R e to three sets of of experimental values: the three 
Peterson et al. |[25l p w values, about 70 Av=l H-W ratios from 
their own experimental spectrum, and Av=2 to Av=l ratios for 
nine upper J levels. The uncertainties of their Av= 1 H-W ra¬ 
tios were larger than for our Av=2 H-W ratios, possibly partly 
due to the Av=2 region being much less contaminated by other 
lines. For the lower vibrational bands, our calculated Av= 1 H- 
W ratios match their observed values about as well as theirs do. 
These differences increase for the higher vibrational bands, but 
the overall weighted root mean square error using our calcu¬ 
lated Einstein A values compared to those of HITRAN is only 
1.41 times higher (for Fn), and 1.22 (for F 22 ). For comparison, 
for the Av=2 H-W ratios, these errors are 3.27 and 4.02 times 
higher using the HITRAN Einstein A values compared to using 
ours. 

A figure equivalent to Figure [3] was produced using the final 
data (Figure|6|, which shows a good match to the observed data. 
Our calculated lifetime (59.13 ms) is also much closer to the 
experimental measurement of van de Meerakker et al. l60l than 
that calculated using the HITRAN DMF (56.6 ms; Table |4j. 


—■— HITRAN 

Solid lines - F1-F1 

—*— Our calculations 

Dotted lines - F2-F2 

—<•— Observed 




Figure 6: Observed and calculated H-W effect in the OH X 2 n, (2,0) band, using 
the final DMF and calculation method. See Figure[3]for full description. 


6. Final Line List, Lifetimes, and Band Strengths 


A line list was produced for all possible rovibrational and 
rotational transitions within the X 2 n state, for v up to 13, and J 
up to between 9.5 and 59.5, depending on the band. Transitions 
are reported between J levels that are a few higher than those 
observed for each vibrational level. Vibrational lifetimes were 
calculated for all of the v levels observed, and are shown in 
Table [5] 

Einstein A v - v values have been calculated for all reported vi¬ 
brational bands, and these are available in the online supple¬ 
mentary material. They are calculated by summing over the 
Einstein Aj'j values for all possible transitions from J' = 1.5, 
F'= 1, for each band. The Einstein A V ' V values have also 
been converted into vibrational band oscillator strengths (/ V ' V 


Table 5: Lifetimes of the first 13 vibrational levels of the OH X 2 n state. 


V 

Lifetime (ms) 

0 

59.13 

l 

31.20 

2 

21.40 

3 

16.07 

4 

12.54 

5 

10.01 

6 

8.12 

7 

6.75 

8 

5.77 

9 

5.11 

10 

4.76 

11 

4.73 

12 

5.16 
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-values) using the equation |f83l 


1 (2 - 6 0 A') 

/ v , v = 1.4991938 ^ /' ’ a v , v , (20) 

v z (2 - 6 0 ,a") 

where A' = A" = 1. 

7. Partition Function 

The partition sums currently provided with the HITRAN 
database are those from TIPS2011 code described by Laraia 
et al. If84l . The values are given between 70 and 3000 K. Con¬ 
sidering that OH is observed in both very cold and very hot 
environments, we decided to extend this range to 5-6000 K. 
As we now have an extended number of calculated energy lev¬ 
els, we calculated the partition function using direct summation 
over the levels derived in Section [2] Our values agree well with 
those calculated using TIPS2011 and those provided in the JPL 
spectroscopic catalogue (based on ref. Il85l ). 

The hydrogen 1=1/2 hyperfine structure is taken into account 
simply by including a twofold degeneracy for each level, con¬ 
sistent with the approach taken by HITRAN. As the hyperfine 
splitting is very small, this makes a negligible difference to the 
calculated values. Figure [ 7 ] shows the relative differences be¬ 
tween partition sums calculated here, and those calculated using 
TIPS2011 and through direct summation of the energy levels 
given in the JPL database (which include hyperfine splitting for 
every level). In general there is good agreement; the discrep¬ 
ancy at low J is an artefact of the interpolation scheme used in 
TIPS2011, which is based on values tabulated at 25 K intervals. 
The growing divergence from the JPL values at higher tempera¬ 
tures appears because the JPL catalogue includes only first three 
vibrational levels. However, in this work we now provide the 
data between 5-6000 K, which is a significant and important 
extension to the previously available partition sums. The list of 
partition sums is given in the supplementary material. 



Figure 7: Relative differences between the partition sums calculated in this 
work, and those calculated using TIPS2011 and through direct summation of 
the energy levels given in the JPL database. TW stands for ’’this work”. 

8. OH and the Oxygen Abundances of the Sun and Stars 

8.1. General Remarks 

OH rovibrational bands can be important features for deter¬ 
mination of oxygen abundances in cool stars (those with effec¬ 
tive temperatures T e ff < 5500 K). They can be attractive al¬ 
ternatives to the visible-region 6300, 6363 A [O I] red lines, 
which are often weak and blended with other atomic/molecular 
transitions, and to the 7770 A high-excitation O I triplet lines, 
whose strengths are subject to significant departures from lo¬ 
cal thermodynamic equilibrium (LTE). Here we report prelim¬ 
inary oxygen abundance analyses from OH lines for a handful 
of stars, and compare our results to published abundances in the 
literature. 

We applied the new OH rovibrational line list to high- 
resolution spectra of the Sun and three cool giant stars in the 
1.5-1 .8 yum wavelength region (the astronomical photometric H 
band). Many lines of the Av=2 system occur in the H band, 
but a large fraction of them are blended with transitions of neu¬ 
tral atomic species, and especially with molecular CN and CO 
bands. Therefore for the present exploratory analysis we lim¬ 
ited our line list in this spectral region to 15 relatively clean 
transitions used in a pioneering study of the Av=2 band system 
( 86 ). 

Several large telescopes now have instruments capable of 
obtaining high resolution (R = A/AA - 20,000-80,000), 
high signal-to-noise ( S / N ) H -band spectra, e.g. VLT CRIRES 
1571 . SDSS APOGEE 1551 . and newly-commissioned McDon¬ 
ald IGRINS 11591 . These facilities will yield a rich set of cool 
stars with easily detectable OH lines. The four stars chosen for 
analysis here represent different temperatures, gravities, metal- 
licities and Galactic population memberships, and they have 
published oxygen abundances. These stars offer differing chal¬ 
lenges for OH analyses. In Figure [ 8 ] we display a small spectral 
interval in each of our program stars. This wavelength range 
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Table 7: Parameters and O Abundances from Individual OH Lines 


x 

E 

Q) 

> 


03 

Q> 


K 




16885 16890 16895 16900 16905 16910 


A 

/jm 

X 

eV 

log (gf) 

Sun 3 
km s -1 

Arcturus 

M67 b 

M71 c 

1.52785 

0.205 

-5.435 


8.72 

8.76 

8.46 

1.54091 

0.255 

-5.435 

8.74 

8.72 

8.76 

8.44 

1.55687 

0.299 

-5.337 

8.74 

8.67 

8.78 

8.45 

1.60527 

0.639 

-4.976 

8.79 

8.69 

8.82 

8.44 

1.61921 

0.688 

-4.957 

8.79 

8.69 

8.75 

8.44 

1.63681 

0.730 

-4.858 

8.74 

8.65 


8.48 

1.64560 

0.609 

-5.108 

8.69 

8.68 



1.65345 

0.781 

-4.806 

8.59 

8.67 


8.42 

1.66054 

0.981 

-4.816 

8.69 

8.67 

8.77 

8.42 

1.66559 

0.681 

-5.069 


8.70 

8.83 

8.47 

1.68722 

0.759 

-5.032 

8.69 


8.78 

8.42 

1.68862 

1.059 

-4.720 

8.74 

8.62 

8.75 

8.41 

1.69042 

0.896 

-4.712 


8.66 

8.75 

8.45 

1.69092 

0.897 

-4.712 


8.69 

8.80 

8.45 

1.76189 d 

1.243 

-4.454 






mean 

8.72 

8.68 

8.78 

8.44 



+/- 

0.02 

0.01 

0.01 

0.01 


sigma 

0.06 

0.03 

0.03 

0.02 


#lines 

10 

13 

11 

13 


a all abundances are in units of log e 
b star 2M08490674+1129529 in open cluster M67 
c star 2M19533986+1843530 in globular cluster M71 
d 1.76189 /im from ref. [86] was unable to be used in any of our program 
stars, but is listed here for completeness 


Wavelength (A) of elements A and B with respect to their solar ratio is written 

as [A/B] = log (Na/Nb)* - log (N A /N H ) Q ; and (c), metallicity 
will be taken to be the [Fe/H] value. 


Figure 8: Observed spectra of the four program stars in the wavelength interval 
covered by Figure 1 in ref. (Hi. Vertical dotted lines denote wavelengths of 
the OH lines. Telluric lines have been removed from the stellar spectra (panels 
b-e) but not from the solar photospheric spectrum (panel a). 

was chosen to match that in Figure 1 of Melendez 1861 . The 
OH lines are strong in Arcturus and the cool giants of M67 and 
M71, but are challenging to detect in the other stars. 

Lines of the OH Av=l system have been studied as part of 
a comprehensive investigation of the solar oxygen abundance 
l90l . The Av=l lines are much stronger than the Av=2 ones, 
but they occur in the 3-4 pm spectral domain (the photometric 
L band). This bandpass is difficult to access with ground-based 
high-resolution spectroscopy because the thermal background 
is higher and the telluric HiO absorption is more time-variable 
in the L band than in the H band. It is challenging to construct 
efficient high-resolution spectrographs for the L band. There¬ 
fore relatively few chemical composition studies have featured 
L-band data. However, high-quality spectra are available for 
the Sun and the mildly metal-poor giant Arcturus, and we have 
applied the new line lists to Av=l lines in these two stars. 

In jj |8.2| we describe the analyses of each star in turn. Standard 
stellar spectroscopic notation are employed here: (a) the “abso¬ 
lute” abundances of a given element A in a star is defined as log 
e(A) = log ((Va/(Vh) + 12.0 ; (b) the relative abundance ratio 


8.2. Abundances in Individual Stars 
8.2.1. The Solar Photosphere 

We employed the 15 lines given in Table 1 of Melendez 
lf86l but did not adopt their equivalent widths (EW) for our 
solar OH analysis. These lines are all extremely weak: 
0.7 mA < EW < 1.4 mA, with (EW) « 1.1 mA= l.lxlO -7 pm. 
This translates to a mean reduced width of ( RW ) = (EW)/A 
^ -7.2. Such lines are at the limit of reliable EW measurement, 
and given the potential for contamination by other solar photo- 
spheric and telluric absorptions, we opted to derive an oxygen 
abundance for each line by matching synthetic and observed 
spectra. Our analyses followed closely those of previous pa¬ 
pers in this series on MgH ll9Tfl , C 2 If92] , and CN (93]; here we 
briefly summarize those methods. 

We synthesized a small spectral range, usually 12 A, sur¬ 
rounding each OH line, using the LTE stellar spectrum analysis 
code MOOG (940 Assembly of the synthesis transition list 
began with the atomic line database of Kurucz |95], adding or 
updating the line parameters of atomic species that have had re¬ 
cent lab analyses by the Wisconsin atomic physics group ( l96l 
and references therein for Fe-group elements; 197] and refer¬ 
ences therein for neutron-capture elements). Molecular transi¬ 
tions of CO l95l , CN 1931 and OH [this study] were merged 


3 Available at http://www.as.utexas.edu/~chris/moog.html 
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Table 6: Stellar Model Parameters and Abundance Summary 


Star 

Te// 

K 

log(g) 

Vmicro 

km s 1 

[Fe/H] 

log E(C) 
adopted 

log e(N) 
adopted 

log e(O) 
[O 1] 

log 6(0) 
OH 

[O/Fe] 
[O I] 

[O/Fe] 

OH 

Sun 

5780 

4.44 

0.85 

0.00 

8.43 

7.83 

8.69 

8.72 

+0.00 

+0.03 

Arcturus 

4286 

1.66 

2.00 

-0.52 

8.02 

7.66 

8.62 

8.68 

+0.45 

+0.51 

M67 a 

4623 

2.44 

1.62 

+0.06 

8.38 

7.94 


8.78 


+0.03 

M71 b 

4314 

1.48 

2.00 

-0.77 

(7.36) 

(7.06) 


8.44 


+0.52 


a star 2M08490674+1129529 in open cluster M67 
b star 2M19533986+1843530 in globular cluster M71 


with the atomic data. Since these line lists were also to be ap¬ 
plied to the spectra of giant evolved stars, we include the major 
CNO isotopologues: 12 CN, 13 CN, 12 CO, and 13 CO, 

For the solar analysis we adopted the photospheric abun¬ 
dances recommended by Asplund et al. {98], in particular 
log e(C) = 8.43, log e(N) = 7.83, log e(O) = 8.69, and 
12 C/ 13 C = 89. To maintain consistency with previous papers 
in this series, we used the empirical model photosphere of Hol- 
weger and Muller |99] . Model atmosphere parameters for the 
Sun and all other stars are given in Table [6] The observed spec¬ 
trum was that obtained by by L. Delbouille, G. Brault, J.W. 
Brault, and L. Testerman at Kitt Peak National Observatory^] 
With the model, line list and observed spectrum inputs we cre¬ 
ated synthetic spectra and altered the oxygen abundance until 
best synthesis/observation matches were obtained. Results for 
each line and the mean abundance statistics are listed in Table[7] 

Our OH solar photospheric oxygen abundance, log e(O) = 
8.72, is only 0.03 dex larger than that recommended by As¬ 
plund et al. |[98] from multiple oxygen-containing species (sig¬ 
nificantly weighted by the [O I] 6300 A transition). Given the 
extreme weakness of the OH Av=2 solar lines, we regard the 
offset as negligible. We also examined the photospheric /.-hand 
spectrum, finding a few very weak and unblended OH Av= 1 
transitions. From these we estimate (log(e(0)) ^ 8.75, in rea¬ 
sonable agreement with the results from the Av=2 transitions. 

8.2.2. Thick Disk Red Giant Arcturus 

The very bright mildly metal-poor star Arcturus has been 
studied many times at high spectral resolution. Our recent study 
of the CN red and violet systems |93l featured Arcturus; that 
paper gives references to previous analyses. Its OH lines are 
strong (panel (b) of Figure [8] and straightforward to analyze. 
We used model atmosphere parameters derived by Ram et al. 
l92| to generate synthetic spectra, and compared them to the 
Arcturus spectral atlas II1001 . The resulting mean O abundance 
(Table [6] of log e - 8.68 translates to a relative overabundance 
of [O/Fe] = +0.51. This value is slightly larger than that derived 
from the [O I] lines by Sneden et al. {93], but is in excellent ac¬ 
cord with that derived by Ram et al. (92). 

8.2.3. Red Giants in Galactic Open and Globular Clusters 

Members of star clusters are tempting targets for chemical 

composition studies since their atmospheric parameters T e ff 


4 Available at http://bass2000.obspm.fr/solar_spect.php 


and log g can often be reliably estimated simply by broad-band 
photometry. APOGEE^ 11011 is a dedicated high-resolution 
spectroscopic survey of more than 100,000 stars in the H band. 
The APOGEE instrument covers most of the H band at a resolv¬ 
ing power of R — 20,000, and currently is producing publicly- 
available spectra for both field and cluster red giant targets. 
From a list kindly supplied by Gail Zasowsky, we acquired 
APOGEE spectra of solar-metallicity old open cluster M67 star 
2M08490674+1129529, and mildly metal-poor globular clus¬ 
ter M71 star 2M19533986+1843530. The APOGEE database^] 
gives derived atmospheric quantities for these stars that we 
adopt here (Table [6]. These stars do not have previously pub¬ 
lished O abundances, but the values derived here are consistent 
with expectations for their clusters: (a) [O/Fe] - 0.0 for the the 
metal-rich M67, and (b) [O/Fe] - +0.5 for the lower-metallicity 
M71. Detailed comparison of O abundances from OH and [OI] 
features should be done in the future, along with derivation of 
O abundances in all APOGEE targets in as many clusters as 
possible. 

In summary, the new OH line lists reported in this paper yield 
O abundances that are in excellent agreement with those deter¬ 
mined from visible spectral region features. The Einstein A 
values used in these calculations are around 13% lower than in 
HITRAN, which translates to a log(gf) change of -0.06 (from 
HITRAN to this work). If the calculations were performed with 
the HITRAN values, the resulting abundances would therefore 
be approximately 0.06 dex lower. For the Sun, this would result 
in log e(O) = 8.66, in equally good agreement with the rec¬ 
ommended values as our result. The HITRAN value for Arc¬ 
turus would be log e(O) = 8.66, equal to the value shown in 
Table [6] (calculated by Sneden et al. (93]), but our value is in 
better agreement with another recently calculated value (8.76 
+0.17; reported by Ramirez and Allende Prieto Il02l ). Pre¬ 
liminary analyses of OH lines appearing in new high-resolution 
IGRINS {89] spectra of several evolved stars have also shown 
better matches to previous abundances for our values than those 
of HITRAN. 

9. Summary and Conclusion 

New absolute line intensities have been calculated for the OH 
Meinel system that we believe to be the most accurate avail- 


5 APOGEE: The ApacheePoint Observatory Galactic Evolution Experi- 
ment; see https://www.sdss3.org/surveys/apogee.php 

6 http://data.sdss3.org/bulkIRSpectra/apogeeID 
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able. The absolute line intensities in common use are based 
on a DMF from 1990 l26l . and using the newly calculated 
UCCSD(T)/ACPF mixed DMF, a better match to an experimen¬ 
tal lifetime l60l and Av=2 H-W ratios is observed. Similarly 
good agreement is seen with the experimental // v values, but 
slightly worse agreement with the Av= 1 H-W ratios observed 
by Nelson et al. l26l . The new line list includes transitions for 
v up to 13, and J up to between 9.5 and 59.5, depending on the 
band, and is available in the online supplementary material. 

The “transformation equation”, that converts transition MEs 
from Hund’s case (b) to (a) (and between LEVEL and PGOPHER) 
was changed to include AX^O MEs. For molecules with a small 
H-W effect, the effect of this change is very small. It is impor¬ 
tant when the H-W effect is large such as for OH, but also for 
NH, for which we recently published a list of intensities 1641 . 
The NH list will be updated with the use of this adjusted method 
in the very near future. 

Partition sums have been calculated using the energy levels 
calculated in this work, and are available in the online supple¬ 
mentary material. They cover the temperature range 5-6000 K, 
which is an increased range compared to what was previously 
available in HITRAN (70-3000 K). This will enable the calcula¬ 
tion of line intensities in cold environments such as interstellar 
clouds, and hot environments such as stars. 

The absolute intensities for 15 Av=2 transitions have been 
used to calculate O abundances in the Sun and three other stars, 
which are compared to abundances recommended in the litera¬ 
ture where available. The OH-based O abundances for the Sun 
and Arcturus are in agreement with the values derived from the 
[OI] 6300 A lines, and those for the two cluster targets are con¬ 
sistent with their expected O abundances. Additionally, our pre¬ 
liminary analyses of OH lines appearing in new high-resolution 
IGRINS ||89l spectra of several evolved stars, including an ex¬ 
tremely metal-deficient star, yield very reliable O abundances. 

The visible spectral region [O I] lines can be difficult to 
work with in many stars, due to contamination by other 
atomic/molecular stellar features and by telluric O 2 absorption 
and night-sky [O I] emission. Given good high-resolution H- 
band spectra, the OH lines in some cases will easily be able to 
serve as the primary O abundance indicators, particularly for 
target stars that are heavily reddened. 

The line list produced in this work will be useful in the fields 
of atmospheric science, astronomy, and combustion science. 
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